On the mobility of extended bodies in viscous films and membranes 
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We develop general methods to calculate the mobilities of extended bodies in (or associated with) 
membranes and films. We demonstrate a striking difference between in-plane motion of rod-like 
inclusions and the corresponding case of bulk (three-dimensional) fluids: for rotations and motion 
perpendicular to the rod axis, we find purely local drag, in which the drag coefficient is purely 
algebraic in the rod dimensions. These results, as well as the calculational methods are applicable to 
such problems as the diffusion of objects in or associated with Langmuir films and lipid membranes. 
The methods can also be simply extended to treat viscoelastic systems. 



I. INTRODUCTION 

The motion of objects in biomembranes is important 
in many cellular processes. These objects are in many 
cases extended, macromolecular inclusions such as pro- 
teins 0, Q or "rafts" [1 El of lipids. Thus, these can 
often be viewed as macroscopic objects moving in a con- 
tinuum fluid environment. Studies of the motion of such 
inclusions in amphaphilic films |gj, |^ and cell membranes 
have a long and interesting history. There are discrep- 
ancies in the early literature on protein diffusion in cell 
membranes, for instance, because of confusion over the 
applicability of three- versus two-dimensional diffusion 
to this case 0, E|- One of the most important contri- 
butions in this area was that of Saffman [8|, who noted 
that the two-dimensional motion in real membranes in- 
duces flows in the surrounding bulk (three-dimensional) 
fluid. He showed that the linearized Stokes-law does not 
describe the motion of inclusions in or bound to mem- 
branes. Rather, this represents a singular perturbation 
problem, and the drag coefficient is not a linear function 
of its size and the viscosity. In fact, the dependence is 
logarithmic, and thus the mobilities and diffusion coeffi- 
cients of proteins in membranes are nearly independent 
of the object's size, in practice. 

Furthermore, there is a length-scale determined by 
the ratio of membrane and fluid viscosities that deter- 
mines the degree to which the dissipation is predomi- 
nantly two- or three-dimensional 0>|ll>l3>G3- We call this 
length £q = 7] m /r)i, where r] mt { are, respectively, the mem- 
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brane (two-dimensional) and fluid (three-dimensional) 
viscosities. A similar length, constructed from the two- 
dimensional shear modulus and the fluid viscosity deter- 
mines the scale of deformations controlled by in-plane 
versus fluid stresses in the case of elastic or viscoelastic 
membranes ^3 E El- Here, we consider the motion 
of extended objects of large aspect ratio in such quasi- 
two-dimensional systems of, e.g., a viscous or viscoelas- 
tic membrane surrounded by a viscous solvent. In such 
cases, both short- and large-scale dissipation can play a 
role. We examine in detail the motion of rod-like inclu- 
sions, but we present a general scheme for the calculation 
of the translational and rotational mobilities of arbitrary 
extended bodies. 

Even for the case of the motion of simple, rod-like ob- 
jects in a bulk Newtonian fluid (i.e., in three dimensions), 
the situation is subtle fA ITiL Hq . For instance, for either 
infinitely long rods in three dimensions or for the motion 
of point-like objects in two dimensions, there is no true 
low Reynolds number regime (i.e., linear hydrodynam- 
ics) [Batchelor 1967]. Specifically, there is no linear drag 
coefficient for such a system: the drag force does not de- 
pend linearly with velocity. For finite length rods and 
small enough velocities, however, there is a linear drag 
coefficient. The drag coefficients for motion parallel and 
perpendicular to the rod axis are given by 

Ci = 2C|| = 4tt?7/ ln( AL/a) 

per unit length, where L is the rod length, a is its radius, 
and A is a constant of order unity, depending on the 
precise geometry. 

We examine here a variant of this hydrodynamic prob- 
lem, in which the rigid rod lies in a two-dimensional in- 
terface that is viscously coupled to a bulk, Newtonian 
liquid phase. The generalization of this problem to the 
motion of such a rod embedded in a viscoelastic film is 
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straight-forward. There are many physical realizations 
of motion in viscous and viscoelastic films. These in- 
clude, for instance, the motion of extended, membrane- 
associated proteins in or on the surface of lipid bilay- 
ers p], Q- Our work here is also motivated by recent 
experiments that have demonstrated the possibility of 
constructing and making quantitative measurements on 
viscoelastic films that closely resemble cellular structures 
such as the actin cortex Driven motion of rods in 
viscous/viscoelastic films has also been used to determine 
rheological properties, e.g., of monolayers [l^. In addi- 
tion, the understanding of this problem will allow the 
quantitative interpretation of the thermally excited an- 
gular fluctuations of microscopic rod-like particles (such 
as fd virus) embedded in such systems. Such hydrody- 
namic studies may also shed light on the dynamics of 
inclusions and transmembrane protein complexes in fluid 
cell membranes. Furthermore, we note that calculational 
methods developed in this work allow one to compute 
the hydrodynamic drag of irregularly shaped objects em- 
bedded in the film. Such objects can be e.g. fractal 
aggregates [l6|. or lipid rafts in call membranes. 

We find two principal results: (i) for small objects 
(specifically, for which L <g; £ ), the drag coefficients be- 
come independent of both the rod orientation and aspect 
ratio; and (ii) for larger rods of large aspect ratio, Q±_ 
becomes purely linear in the rod length L — i.e., the drag 
becomes purely local. In contrast, we find that the well- 
established three-dimensional result C|| = 2^77/ ln(AL/a) 
applies for parallel motion in the film, provided that 
L > 4- Here, however, the effective rod radius be- 
comes £q rather than the physical radius a, when a <C £o- 
Closely related to (ii), we find that the rotational drag 
(equivalently diffusion constant) depends purely alge- 
braically on the rod length. 

We consider only the most simple class of rod geome- 
tries in that we assume the rod to have a circular cross 
section in the plane perpendicular to its axis. We take 
the radius of that circle to be a while the length of the rod 
is given by L. The geometry of this rod of length L and 
radius a is then parameterized by only one dimensionlcss 
ratio, the aspect ratio p defined by 



9 



L 

2a' 



(1) 



The long axis of the rod may be assumed to be termi- 
nated by spherical end caps although, as will be seen be- 
low, our calculation is not sensitive to the fine structure 
of the ends due to our introduction of a short distance 
cutoff in the problem. We will argue below, however, 
that for the case of large aspect ratio rods, the detailed 
structure of the end caps will have a negligible effect on 
the overall drag. 

The fundamental distinction between the well-known 
result for the hydrodynamic drag on a rod in a homoge- 
neous, three-dimensional viscous material and our result 
for the drag on a rod embedded in a viscous membrane 
coupled to a viscous, three-dimensional fluid is the ap- 



pearance of the length scale £q. This length enters since 
the ratio of a two-dimensional interfacial viscosity rj m 
and the usual viscosity 7/f of the bulk fluid has dimen- 
sions of length, 



4 = 



(2) 



The denominator, J^rj is the sum of the viscosities of 
the fluid above the interface (the superphase) and fluid 
below the membrane (the subphase), in the general case 
of two different viscosities (or one vanishing one, in the 
case of Langmuir films). For a membrane embedded in a 
uniform fluid of viscosity 77, this is ^2 V = 2f?- In general, 
however, we may have two distinct nonzero viscosities. 
For a Langmuir monolayer we have ^2 V = ?7j since we 
can neglect the viscosity of air. From here on, however, 
we shall simply use the length £q and rj^ for the sum of the 
viscosities of the two bulk fluids surrounding the mem- 
brane. Throughout this paper we measure all lengths in 
terms of the fundamental length £q, except where explic- 
itly stated to the contrary. 
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FIG. 1: The upper figure shows the film from the top looking 
down. The rod is shown as the black line and the film is 
shown in blue. The rod is assumed to be embedded in the 
film, but not in the bulk subphase below it. This subphase is 
shown in the lower figure, which pictures the system from a 
side view. In this paper the Newtonian subphase is assumed 
to be infinitely deep. 



There are three independent drag coefficients or parti- 
cle mobilities to determine in the problem. The mobility 
tensor is the inverse of the drag coefficient tensor. For 
translational motion, the mobility tensor /Xy of the rod 
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in the membrane is defined by: 
„.rod 



HijF. 



j ■■ 



(3) 



where vj 0( ^ is the i th 



rod and Fa is the j 



th 



component of the velocity of the 
component of the force applied to 



the rod. By in-plane rotational symmetry combined with 
the n — ► — h symmetry of the rod, the mobility tensor 
must take the form: 



M-L {h. 



(4) 



Here fj,± and fj,u are the mobility of the rod dragged per- 
pendicular to and parallel to its long axis respectively. In 
addition to these two independent translational mobili- 
ties, there is also one rotational mobility, fx TOt linking the 
angular velocity of the rod to the torque applied to that 
rod about its center of inversion symmetry. 

It should be noted that there is no hydrodynamic 
torque acting on the rod when it is dragged by any force 
acting through the center of inversion symmetry of the 
rod. We refer to these forces as "central." This can be 
seen from the following argument: Torque in the two- 
dimensional plane is a pseudoscalar that must be odd 
under time reversal. Such a pseudoscalar can only be 
built out of the antisymmetric combination of velocity 
vector of the rod v and a vector along the axis of the rod 
n. That combination e a pv a np must also be symmetric 
under n — ► — n since the rod is symmetric upon such 
transformation. Thus the only available pseudoscalar is 
disallowed. Therefore there is no rotational motion gen- 
erated by this class of central forces and the rotational 
degree of freedom decouples from the two translational 
degrees of freedom. We further note that boundary con- 
ditions that break rotational symmetry or the application 
of a force at a point other than the center of symmetry 
violate the above assumptions and thus allow a coupling 
of the rotational and translational motion of the rod. 

We approach the solution of this problem by two com- 
plimentary means. In the first part of the calculation, 
section Til Al we approximate the continuous rod by a se- 
ries of discs, in analogy to the Kirkwood approximation 
used in the calculation of the drag on a rod in a uniform, 
three-dimensional viscous environment [l^. This calcu- 
lation allows one to incorporate the details of the shape of 
the rod in the resulting drag coefficient. Here the shape of 
the rod is parameterized in terms of its dimcnsionless as- 
pect ratio. This method, however, becomes numerically 
difficult in the limit of large aspect ratios, i.e., for very 
long, thin rods. In this limit we may proceed by a sec- 
ond approximation that assumes the rod to be infinitely 
thin — i.e., L>a. Here, we also restrict our attention to 
the limit a <C (q. 

The outline of the remaining parts of this paper are 
as follows: In section III Al we develop the first of two 
calculational methods for determining the drag on rod- 
like objects embedded in the membrane. Then in section 
IIIBI we demonstrate a second approach to determine the 
drag on a rod. This approach is optimized to work in 



the limit of high aspect ratio rods and compliments the 
first approach which is most efficient for rods of small 
aspect ratio, i.e., less elongated objects. The reader who 
wishes to examine the results of these calculations with- 
out delving into their details may skip to section lTTll whcrc 
the drag coefficients for translational motion parallel and 
perpendicular to the long axis of the rod are computed 
for a variety rod geometries. In addition we explore the 
rotational drag on these rods there. Finally, we conclude 
and discuss our results in section llVl 



II. THE RESPONSE FUNCTION 

A. By the Kirkwood approximation 

To incorporate the correct dynamics of this coupled 
system, we use the results of our previous calculation|l9| 
for the velocity response of the membrane at a position 
x due to the application of a force localized at the origin. 
Specializing to the case of a purely viscous (as opposed 
to a viscoelastic) film, we can write the membrane ve- 
locity response to a point force localized at x' using the 
previously defined response function: 



v a (x) = a af3 (x - x') fp (x') 



(5) 



in a closed form as 



a a p (x) 



(|x|) x a xp + ct|| (|x|) [6 a p - x a xp] , (6) 



where the scalar functions aii , a± of the distance between 
the point of the force application and the measurement 
of the velocity field are given by: 



1 



(7) 
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l Ul (z)-js-~(Y (z)+Y 2 (z)) 



and 



-iuiocx(x, to) = 
1 



(8) 



47T?7„ 



ttHoO?) - -Hi(z) 



(Y (z)-Y 2 (z)) 



where the are Struve functionspCI]. and the Y v are 
Bessel functions of the second kind. We note z = |x| /£o 
is the distance between the point of force application 
and the membrane velocity response measured in the flat 
membrane in units of £q- At this point, it is important 
to distinguish between the two response functions intro- 
duced: the mobility tensor of the rod gives the veloc- 
ity response Vi of the rod given a total force Fj acting on 
it through its center. The response function ct^x — x') 
gives the velocity response Uj(x) of the two-dimensional 
fluid at the position x due the application of a point force 
Fj(x') at another point x' in the two-dimensional fluid. 
The main purpose of this and the next section of the pa- 
per it to derive the former response function for the rod 
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from the latter response function for the fluid, which we 
have previously calculated pjj. 

We note also that, due to the linearity of the under- 
lying low-Reynolds number hydrodynamics, the velocity 
response produced at some point in the membrane by a 
collection of point forces is simply the sum of the response 
functions appropriate to each point force individually: 

JV 

v a (x) = 7^ a a /3 (x - x n ) fp (x w ) (9) 

n=l 

where n indexes the N point forces located respectively 
at locations x„ in the film. Clearly the sum may be 
converted into an integral for case of a continuous force 
distribution; we will examine this limit for the case of an 
infinitely thin rod of finite length in the next section. 

Using this superposition principle we may determine 
the effective drag on a rod by employing the two- 
dimensional analogue of the Kirkwood approximation 
used to calculate the hydrodynamic drag on a rod in three 
dimensions. Specifically, we replace the rod of length L 
and cross-sectional radius a by a set of n + 1 disks of 
radius a and intersphere separation b chosen so that the 
total length of the rod is preserved, i.e. L = nb + 2a. 
See figure El In this way, the aspect ratio of the rod, 
p = L/a can be fixed. The number of disks, of course, 
can be varied, however, we will always choose that num- 
ber to be the maximal number consistent with a given 
aspect ratio and the non-interpenetrability of the disks. 
We have also shown that computed hydrodynamic drag 
on the rod is only weakly dependent on the number of 
disks (or equivalently on the inter-disk separation). 

Our strategy for computing the drag on the rod in- 
volves setting the rod in uniform motion with unit veloc- 
ity by imposing some set of forces f"), i — 1, . . . ,n + 1 
on the n + 1 disks that make up the rod. Clearly the to- 
tal force on the rod, which is equal to the effective drag 
coefficient is simply the sum of those forces: 

n+l 

Fa = fa = ^ ( 10 ) 

i=l 

Using Eq. JSJ we can compute the velocity field for a 
given collection of point forces. However, to determine 
the effective drag on the rod, one needs to insist that all 
the beads have the same given velocity and determine 
the forces applied to them to ensure this result. 

For definiteness, we first discuss the drag on the rod 
when moving perpendicularly to its long axis. We may 
then directly write out the analogous solution for motion 
parallel to the long axis of the rod. Motion of the rod 
along any arbitrary direction in the plane relative to its 
long axis can be obtained from these two results using 
the linearity of problem. Thus these solutions span all 
possible linear motions of the rod. We will return to 
rotational motion shortly. 

To calculate the drag on the rod moving perpendicular 
to its long axis (see figure P), we first apply Eq. ©to 




FIG. 2: Approximating the continuous rod by a series of balls 
of radius a with center-to-center separation b. The number 
of balls and their separation is chosen so that collection of 
balls has the same length L as the original rod. For the best 
approximation to the original rod, we maximize the number 
of balls for a given aspect ratio of the original rod. 



determine the velocity of all n + 1 spheres making up the 
rod as a function of as yet unknown forces. We insist 
only, based on the symmetry of this problem, that these 
forces are also perpendicular the rod's axis. Thus, 

n+l 

^>=$>f /« (11) 

i=l 

where we have suppressed the vectorial indices, defined 
to be the velocity of the sphere at position x^, 
and rewritten the mobility tensor using Eq. JSJ) and the 
definition: 

af =oi(x«-xW). (12) 

The solution for the forces on the individual beads and 
thus, using Eq. (JTDJ , the total force on the rod and equiv- 
alently the hydrodynamic drag is then obtained by in- 
verting the (n+l) x (n+l) matrix of response functions 

a± ■ We call this inverse matrix defined by: 
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The drag coefficient is then 

n+l 

a = £ m™, (i4) 

where we have used Eq. (|1U|I and the fact that each bead 
making up the rod has a velocity of unity. The analogous 
expression from the drag on the rod moving parallel to 
its long axis is then obtained by replacing the terms of 
the mobility matrix, with .Mi since here all the 

bead velocities are parallel to the separation vectors be- 
tween the beads. We then have 

n+l 

C|| = E (15) 

and the complete solution for motion of the rod in the 
plane at an arbitrary angle <fi with respect to its long axis 
(see figure is given by 

C(0) = C|| cos^ + Cisin^. (16) 

It is perhaps not surprising that the drag on the rod 
does not have a simple, closed-form solution in the gen- 
eral case. Here we examine the results of a numerical 
calculation of the these drag coefficients. The results 
depend on two parameters reflecting the relative magni- 
tudes of three lengths involved in the problem. We report 
our results in terms of the two dimensionless numbers in- 
troduced above: the overall length of the rod (i.e., the 
linear dimension of its long axis) measured in units of £q, 
L/£q, and the aspect ratio of the rod L/a. 

B. Thin rod approximation 

As a calculation tool, the above Kirkwood approxi- 
mation has one shortcoming. If one wishes to explore 
the drag on a very thin rod compared to its length,?. e., 
one of high aspect ratio, one is compelled to use a large 
number of disks. In fact, the number of disks is linearly 
proportional to the aspect ratio of the rod being mod- 
elled. Since even numerically inverting a large matrix is 
cumbersome, the penultimate step in the procedure out- 
lined above becomes slow in the high aspect ratio limit. 
Fortunately, there is a complementary approach to the 
Kirkwood scheme that involves inverting a more manage- 
able matrix and is still valid for exploring rods of infinite 
aspect ratio, or specifically, rods of finite length, but van- 
ishing thickness. Since it is precisely in this limit of large 
aspect ratio that the Kirkwood approximation becomes 
intractable, the combination of these two approaches al- 
lows us to study the drag on rod for both large and small 
aspect ratios. 

The thin rod approximation starts from the assump- 
tion that one can take the limit of an infinitely thin rod 
from the start. Thus the velocity field at the point x due 
to a continuous distribution of force densities along the 



rod, f(xx) that lies along the x-axis from x — —L/2 to 
x = L/2 takes the form 

fL/2 

Vi(x) = / oiij (x - px) fj {px) dp, (17) 

J -L/2 

the indices in the above equation represent the usual vec- 
torial indices; there are no disks to count here. There has 
been a simplification introduced in Eq. El The integral 
over the force density, / (x) which in principal extends 
over the entire volume of the rod has been made one- 
dimensional by assuming the rod to be infinitely thin. 
This simplification is possible since the short distance 
behavior of the response function includes only an inte- 
grable singularity. 

After writing Ea.ll7lone is still faced with the problem 
of inverting the equation. After all, the assumed rigid- 
ity of the rod requires that v(x) be constant along the 
rod, but force on any length clement of the infinitely thin 
rod are unknown as they are comprised of both the exter- 
nally applied for and, as yet undetermined internal forces 
of constraint. Our method of solution is fundamentally 
identical to that used in the previous section, we impose 
a unit velocity field on the rod and determine the force 
density required to affect his result, f (x). We use this 
intermediate result by integrating the linear force density 
over the length of the rod to find the total force necessary 
to move the rod with a velocity of unity. This force is 
clearly just the drag coefficient that we sought. 

The inversion of Eq . ITTl Droceeds by first expanding the 
linear force density in Legendre polynomials P n (2x/L). 
From symmetry considerations we can expand the force 
density in even Legendre polynomials, i.e., 

N 

/(a;) = E c 2»- P 2„(2x/L) (18) 

n=0 

for the both the perpendicular and parallel dragging cal- 
culations, while we expand the force density in odd Leg- 
endre functions to study the rotational drag on the ob- 
ject. In what follows, we only describe in detail the case 
of linear drag, although the rotational case is very sim- 
ilar. The coefficients c„ are as yet unknown, since the 
function f(x) is undetermined. In principle since the 
Legendre polynomials form a complete set on the inter- 
val — 1 to +1, any physical force density can be expressed 
as in Eq. 1 181 provided N be taken to infinity. In practice 
we find excellent numerical results even when truncating 
this sum to just the first few terms, taking into account 
the symmetry of the force distribution about the center 
of the rod. We discuss in more detail below the meaning 
of numerical accuracy as used in the present context. 

Using the expansion of the force density from Eq. El 
we can rewrite Eq. El as 

N fL/2 

v(x) = E c 2«/ a(x- zx)P 2n (2z/L)dz. (19) 
n=0 J - L / 2 
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To further simplify the analysis, we relax the requirement 
that the velocity field be equal to unity at all points along 
the rod. Rather we impose this condition at a finite set of 
points < pi < L/2 along the rod. Given that we intend 
to truncate the Legendre function expansion of the force 
density at TV, we can impose the velocity condition at 
a maximum of TV + 1 points without creating an over- 
determined system of equations. Thus we require that: 



v (pix) 



O.fori = 0, ... ,7V 

\ oL/2 



(20) 



1 = 



c 2n / a (xpi - zx) P 2n {^z/L) dz 

J -L/2 



n=0 - - L /2 

for i = 0, . . . , TV. 



(21) 



Now finding the force distribution along the rod requires 



only the inversion of an (TV + 1) x (TV 
whose components are defined by 



1) matrix, A/i, 



L/2 

, a (xpi — zx) P2j(2z/L) dz. (22) 

L/2 



Finally, the total force on the rod is found by recon- 
structing the force density from its Legendre polynomial 
expansion and integrating the resulting expression over 
the entire rod. That force density is given by Eq. ^| 
where the coefficients, are determined by solving the 
set of equations Eqs. 1201 using the inverse of the matrix 
defined in Eq. E2 Thus we find 



JV 



Ck 



(23) 



III. RESULTS 

Figure [3] shows the parallel drag coefficient for a rod 
as a function of its length, L/Iq in reduced units. Two 
finite aspect ratios, as well as the case of infinite aspect 
ratio, i.e., an infinitely thin rod, are shown in this fig- 
ure. As expected based on the above discussion, the 
shape of the rod and thus its aspect ratio does not ef- 
fect the drag on the rod in the limit all dimensions of 
the rod are small compared to l§. Essentially all small 
rod-like particles behave in the membrane as if they were 
infinitely thin. For rods longer than this crossover length, 
the shape of the rod clearly affects its hydrodynamic drag 
as seen by the deviation of the finite-aspect-ratio curves 
from the p — oo curve shown as the solid line in the 
figure, ft should be noted that the finite aspect ratio 
calculations were performed using the Kirkwood approx- 
imation, while the drag on the infinitely thin rod (p = oo) 
is calculated using the thin-rod approximation. 
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Due to the orthonormality of the Legendre polynomials, 
the total force is given entirely by the coefficient of the ze- 
roth Legendre polynomial, cq. Similarly, the total torque 
in the case of rotations is given by the coefficient of the 
first odd Legendre polynomial. 

There remains one more aspect to this calculation. 
How does one determine the accuracy of the approxi- 
mation and how can one optimize that accuracy by se- 
lected the set of points along the rod a;,*, i = 0, . . . , TV at 
which to impose the velocity boundary condition? We 
first judge the accuracy of the result by calculating the 
velocity field along the rod. Ideally the value of the ve- 
locity should be unity. Because we have only imposed 
that boundary condition at a discrete set of points along 
the rod, the velocity field along the rod varies. We find 
that with only eleven points chosen symmetrically about 
the center of the rod (i.e., TV = 5 above), the velocity 
field deviates from unity by at most five percent in most 
cases. However, we find that the drag coefficient (defined 
by the ratio of total force to average velocity) converges 
much more rapidly with TV. Specifically, we find a vari- 
ation of less than a percent in the drag coefficient in all 
cases, so long as TV > 2. Below, we report our results 
for TV = 5. We thus are able to achieve better than one 
percent accuracy with only a 6 x 6 matrix. 



FIG. 3: Parallel drag coefficient of a rod for various aspect 
ratios. 

The analogous calculation can be made for the drag on 
the rod moving in a direction perpendicular to its long 
axis. These results are shown in figure 0] 
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FIG. 4: Perpendicular drag coefficient of a rod for various 
aspect ratios. 

The convergence of the finite aspect ratio rod results 
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towards the infinite aspect ratio result clearly tests the 
consistency of the two numerical approaches to drag cal- 
culation. To study in more detail the aspect ratio depen- 
dence of the drag we plot the drag coefficient for parallel 
motion (figure [5J and perpendicular motion (figure EJ. 
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FIG. 5: The parallel drag coefficient of a finite aspect ratio 
rod as a function of aspect ratio, p. The result for an infinite 
aspect ratio is shown as a horizontal line. The length of the 
rod is 20l • 
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FIG. 6: The perpendicular drag coefficient of a finite aspect 
ratio rod as a function of aspect ratio, p. The result for an 
infinite aspect ratio is shown as a horizontal line. The length 
of the rod is 2Q£ 

In both of these figures [5] and H3 the length of the rod 
was held constant so that L/£q — 20. As discussed above 
particle shape is less relevant for particles with dimen- 
sions less than £q so a rod of significantly longer length 
was chosen to explore the aspect ratio dependence of the 
rod's drag coefficient. To observe the importance of the 
particle size (measured in the natural units of £q) we plot 
the aspect ratio dependence of the drag coefficient of a 
rod of length O.l^o- The drag coefficient of the rod mov- 
ing parallel to its long axis is shown in figure while 
the drag coefficient of the rod moving perpendicular to 
its long axis is shown in figure Q3 

It may be observed from a comparison of the pairs of 
corresponding figures for parallel drag and perpendicular 
drag such asHvs. IBlElvs. © and 03 vs. ^ that the per- 
pendicular drag coefficient is strictly larger than the par- 
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FIG. 7: The parallel drag coefficient of a finite aspect ratio 
rod as a function of aspect ratio, p. The result for an infinite 
aspect ratio is shown as a horizontal line. The length of the 
rod is 0.1£o 
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FIG. 8: The perpendicular drag coefficient of a finite aspect 
ratio rod as a function of aspect ratio, p. The result for an 
infinite aspect ratio is shown as a horizontal line. The length 
of the rod is 0.1£q 



allel drag coefficient for rods of all aspect ratios (greater 
than one, i.e. not disks) and lengths. The magnitude of 
the difference between these two drag coefficients, how- 
ever, depends on the length of the rod compared with 
the natural length, £ . For lengths such that L < £ , the 
two drag coefficients converge to the same value and the 
difference between these coefficients grows monotonically 
with increasing rod length. As an example, the two drag 
coefficients are plotted as a function of rod length for 
infinite aspect ratio rods in figure 

The two coefficients begin to separate at rod lengths 
on the order of £q. This differs not only from the small 
rod limit but also from the case of motion in bulk fluids, 
where the two drag coefficients differ only by a constant 
factor of two. The length dependence here can be under- 
stood by noting that £o sets the natural length scale over 
which the two-dimensional fluid velocity field can vary. 
In detail, what we find is that the parallel drag in the film 
is essentially unchanged from the bulk, three-dimensional 
drag, in that 
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FIG. 9: A comparison of the perpendicular and parallel drag 
coefficients as a function of rod length for a rod of infinite 
aspect ratio. Note that the difference between the two coef- 
ficients is a monotonically increasing function of rod length 
and that the two coefficients begin to diverge at the length 

where the prefactor in the logarithm has been determined 
to within 1%. Comparing with the result for drag of a 
rod in a bulk fluid, the effective radius of the rod is now 
of order £q (for a -C £q). On the one hand, that the 
three-dimensional result is recovered is not surprising, 
given that £q corresponds physically to a length scale be- 
yond which the fluid viscosity dominates the film viscos- 
ity. Furthermore, the corresponding fluid velocity field in 
this case both respects the assumed incompressibility of 
the film, and is the same as that of rod motion in a fluid 
above this length £q. Thus, £o determines the effective 
aspect ratio. The case of perpendicular motion, on the 
other hand, is qualitatively very different. We find 

Ci = 2-nriL. (25) 

Here, the corresponding bulk fluid velocity field in the 
absence of the film is inconsistent with incompressibility 
of the film. Specifically, there is a non-vanishing two- 
dimensional divergence of velocity field restricted to the 
plane of motion for motion perpendicular to the rod axis. 
Hence, although only the bulk fluid viscosity r\ enters this 
expression (to be expected since dissipation is dominated 
at the largest scales by the fluid viscosity), the in-plane 
incompressibility requires that the fluid velocity field ex- 
tends over distances comparable to the largest dimension 
L. This means that the usual hydrodynamic coupling of 
portions of the rod (represented by the logarithm) is not 
present. The result is a drag coefficient purely linear in 
rod length. In other words, the drag is effectively local 
in character. 

Finally, we note the result for the rotational drag coef- 
ficient on the rod which gives the required torque applied 
to rod about the center to generate an angular velocity 
of the rod equal to unity. As discussed in the previous 
sections, this calculation proceeds analogously to those of 
the perpendicular and parallel drag coefficients. We plot 
the rotational drag coefficient divided by L 2 for a rod of 
infinite aspect ratio as a function of the reduced length 
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log (£//„) 

FIG. 10: The solid line is the ratio of the perpendicular drag 
coefficient to the parallel drag coefficient calculated in the 
thin rod approximation. The dotted lines are two different 
asymptotic fits to this curve corresponding to short rods and 
long rods. For short rods this fit is to the constant one con- 
sistent with the Saffrnan-Delbriick result. The large rods fit 
is to a simple logarithm in length as discussed in the text. 
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FIG. 11: The rotational drag coefficient of a rod of infinite 
aspect ratio plotted versus the length of the rod. 

in figure The essential feature of this plot is that 
rotational drag coefficient scales as L 2 for rods smaller 
than £q and then as L 3 for rods longer than this natural 
length. Thus, we find purely algebraic behavior in both 
limits. 



IV. SUMMARY 

Using the response function previous calculated we 
have calculated the hydrodynamic drag on a rod mov- 
ing at low Reynolds number in a viscous film coupled 
to fluid sub- and superphases of arbitrary viscosity. The 
drag coefficient on the rod is a tensorial object with two 
independent parameters that correspond to the drag co- 
efficient of the rod moving along its long axis (parallel) 
and in the direction perpendicular to its long axis. We 
have also computed the rotational drag coefficient. 

These results were calculated numerically using two 
methods with complementary regimes of validity; the 
Kirkwood method, which approximates the rod as a se- 
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ries of non-interpenetrating disks linked together and is 
well suited to calculating the drag coefficients for rods of 
smaller aspect ratio. The aspect ratio is set by choos- 
ing the number of these non-interpenetrating disks to 
make up the rod. For very long, thin rods having higher 
aspect ratios, this method becomes numerically cumber- 
some since it involves inverting an n x n matrix for a 
rod made up of n disks and reaching higher aspect ra- 
tios requires adding more disks. To explore the limit of 
very high aspect ratios, one can perform calculations in 
the infinite aspect ratio limit. These two methods can 
be shown to be consistent numerically; in the limit of a 
large number of disks, we have checked numerically that 
the results of the Kirkwood method approach those of 
the thin rod approximation. 

It is instructive to contrast our results with those for 
the three-dimensional case. In three dimensions, there 
is a length-independent factor of two difference between 
the parallel drag coefficient and the perpendicular drag 
coefficient: 
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The appearance of the logarithm in Eq. signals the 
break down of a purely local drag (or, "free draining" ) . In 
other words, the long-range hydrodynamic interactions 
between various segments of the rod cause the drag on 
the rod to be reduced from simple linear dependence on 
L as would be the case if the hydrodynamic drag on each 
element of the rod were purely local in character and 
thus the total drag additive along the length of the rod. 
Instead, the motion of one part of the rod sets up long- 
ranged fluid flows that effectively drag other parts of the 
rod forward. 

The reduced dimensionality of the flow in the film qual- 
itatively changes this result as show in Eg. 1241 and figure 
|5J From ITU1 it is clear that the two drag coefficients are 
equal in the limit L <c £o (the dotted line for small L /£q 
is simply unity) while in the limit that L ^> £q they differ 
substantially. In fact, as we argue here, we see an appar- 
ent, purely local drag per unit length. Hence, the ratio of 
the two drags for long rods is given simply by the loga- 
rithm described above, as can be seen by the asymptotic 
fit to a logarithm that is illustrated by the dotted line on 
the right of the figure. 

On the one hand, while the dimensions of the rod are 
small (<C £q), the dissipation is governed primarily by the 
film, which is insensitive to orientation and aspect ratio, 
just as it is to size: there is only a weak, logarithmic 
dependence on size in this limit [8(. On the other hand, 
when the rod becomes longer (than £q), the dependence 
on both orientation and aspect ratio becomes stronger. 
Here, the three-dimensional fluid governs the dissipation, 
and we know that orientation and especially size matters 
in this limit. A major difference, however, arises when 
we compare parallel motion with perpendicular motion. 



In the latter case, although the dissipation is governed 
primarily by the fluid, the film (along with its assumed 
incompressibility and no-slip conditions) imposes a very 
different boundary condition on the flow from what we 
would have in a bulk fluid alone. The velocity field v±_ in 
the film satisfies Vi • v± — 0, which is inconsistent with 
the Stokes flow in perpendicular motion. Considered as 
a two-dimensional field, d x v x + d y v y is non-zero in the 
plane of motion in this case. This added condition can 
only increase the drag for perpendicular motion relative 
to that without the film present. In contrast, since this 
boundary condition is consistent with the flow field for 
parallel motion, we expect to find quantitative agreement 
with the parallel mobility in a bulk fluid when the film's 
viscosity becomes irrelevant (small ^o)- 

For the perpendicular motion, not only is the drag in- 
creased relative to that for a bulk fluid, but the depen- 
dence is purely linear, as show in Fig. |5J The linear de- 
pendence is an indication of the absence of the hydrody- 
namic effects described above, which reduce the drag by 
cooperativity of sections along the rod. Here, the drag is 
simply proportional to the length of the rod. This can be 
seen from the additional boundary condition mentioned 
above. Although the dissipation for long rods is governed 
by the fluid viscosity, the characteristic scale for this flow 
is that of the whole rod. Unlike the case of perpendicular 
motion in a simple fluid, where there is a short path of 
order the rod diameter a around the rod, the in-plane 
incompressibility forces the flows to go around the long 
way. Hence, the total absence of the logarithm, and the 
simple, purely local drag proportional to length. This can 
be seen in Fig. I1UI where we show the ratio of the drag 
coefficients is just given by the logarithmic term coming 
from the hydrodynamics of rods in ordinary fluids. Fi- 
nally, we note that these observations also explain why 
the rotational drag is purely algebraic for long rods, since 
rotations exhibit a purely local drag or rod segments per- 
pendicular to the motion. Specifically, as shown in Fig. 
II 11 we find 



(r ~ 0.167T77L J 



(28) 



In summary we have developed a highly adaptable 
framework to compute the drag on irregularly shaped ob- 
jects embedded in a viscous membrane or interface. As 
a demonstration of this method we have computed the 
drag on a rigid rod in this two-dimensional fluid system 
viscously coupled to a fluid subphase and compared our 
results to both the well-known results for the drag on a 
rod in a three-dimensional, viscous fluid, and the result 
for the drag on a disk embedded in a membrane (due to 
Saffman and Delbriick), which is applicable to e.g. the 
diffusion of small transmembrane proteins. Firstly, we 
find, in accordance with the results of previous investiga- 
tors @, El E3 that there is an inherent length scale £ in 
the system set by the ratio of the two-dimensional viscos- 
ity of the membrane to the three dimensional viscosity of 
the subphase fluid. The existence of such a length scale is 
made obvious by dimensional analysis. The importance 
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of this length scale on the drag tensor associated with 
various objects embedded in the membrane has been dis- 
cussed in this work. In brief, for objects with character- 
istic dimensions less than £, the Saffman-Delbruck result 
is recovered from our more general computation. The 
drag coefficient tensor for all such small objects is sim- 
ply that of a small disk. It is isotropic and independent 
of particle size except for logarithmically small correc- 
tions. For objects significantly larger than £, our analy- 
sis shows that drag tensor becomes both size dependent 
and anisotropic for rod-like objects. Based on our cal- 
culations restricted to rod-like objects, it is nevertheless 
clear that these general statements apply to objects of ar- 
bitrary shape as well. In particular for these rods, we find 
that the parallel drag coefficient matches that of a rod 
dragged parallel to its long axis in a three-dimensional, 
viscous fluid. For the same rod dragged perpendicular to 
its length a qualitatively different result is found: the ef- 
fective drag coefficient is logarithmically enhanced versus 
its three-dimensional counter-part due to the breakdown 
of long-range hydrodynamic interactions in the interface 
due to momentum transfer to the bulk, subphase. 

While the heuristic importance of this new length scale 
is clear, we have not yet estimated its size for typical 
membranes. As a starting point, we note that if the mem- 
brane/interfacial viscous were equal to that of the bulk, 
subphase, this length would naturally be the thickness of 
the membrane, i.e. a molecular length. For a short chain 
surfactant monolayer or lipid bilayer this length would 
then be on the order of 1 — 2nm respectively. However, 
it is expected that the internal viscosity of the mem- 
brane/interface is typically much larger than that of the 
(typically aqueous) subphase and this length £ is con- 
sequently multiplied by a factor equal to this viscosity 
enhancement of the membrane/interfacial material. It is 
not unreasonable to suppose that £ ~ 10 — lOOnm. Thus 
we expect to see significant deviations from the Saffman- 
Delbruck result both for lipid rafts and protein aggre- 
gates that surpass such lengths. At the same time, it is 



clear that the standard Saffman-Delbruck result should 
explain the observed mobility of individual transmem- 
brane proteins. 

Further experimental tests of the above theory require 
the analysis of tracer particle diffusion data for a mem- 
brane/interface bound objects of various sizes. Based 
on these calculations, the observation of anisotropic dif- 
fusion constants for these longer rod-like objects would 
be a clear indication of phenomena unexplainable by the 
Saffman-Delbruck analysis. After testing these basic mo- 
bility calculations, one could then use these results to 
do both standard translational microrhcology on mem- 
branes and interfaces ^ij as well as novel experiments on 
rotational microrheologvplj . Such studies will be par- 
ticularly interesting in studying the properties of bro- 
ken rotational phases of lipid monolayers such as hexatic 
phases 22]. In addition one should be able to extend the 
present calculations to explore the implications of mem- 
brane hydrodynamics upon diffusion limited aggregation 
in membranes. Such aggregation processes play an in- 
teresting role in the formation of transmembrane protein 
aggregates, lipid rafts, and colloidal aggregates on large, 
unilamellar vesicles 231. 
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